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Abstract 

We introduce a lattice model for a static and isotropic system of relativistic fer- 
mions. An action principle is formulated, which describes a particle-particle interac- 
tion of all fermions. The model is designed specifically for a numerical analysis of the 
nonlinear interaction, which is expected to lead to the formation of a Dirac sea struc- 
ture. We discuss basic properties of the system. It is proved that the minimum of the 
variational principle is attained. First numerical results reveal an effect of spontaneous 
symmetry breaking. 

1 Introduction 

It is generally believed that the concept of a space-time continuum (like Minkowski space 
or a Lorentzian manifold) should be modified for distances as small as the Planck length. 
The principle of the fermionic projector [1] proposes a mathematical framework for physics 
on the Planck scale in which space-time is discrete. The physical equations are formulated 
via a variational principle for fermionic wave functions defined on a finite set of space-time 
points, without referring to notions like space, time or causality. The idea is that these 
additional structures, which are of course essential for the description of nature, arise as 
a consequence of the nonlinear interaction of the fermions as described by the variational 
principle. More specifically, it was proved that the original permutation symmetry of 
the space-time points is spontaneously broken by the fermionic wave functions [3]. This 
means that the fermions will induce non-trivial relations between the space-time points. 
In particular, one can introduce the notion of a "discrete causal structure" (see the short 
review article [5]). The conjecture is that for systems involving many space-time points 
and many particles, the fermions will group to a "discrete Dirac sea structure", which in a 
suitable limit where the number of particles and space-time points tends to infinity, should 
go over to the well-known Dirac sea structure in the continuum. Then the "discrete causal 
structure" will also go over to the usual causal structure of Minkowski space pQ . 

Hints that the above conjecture is true have been obtained coming from the continuum 
theory. First, our variational principle has a well defined continuum limit [U Chapter 4], 
and we get promising results for the resulting effective continuum theory [TJ Chapters 6- 
8]. Furthermore, rewriting certain composite expressions ad hoc as distributions in the 
continuum, one finds that Dirac sea configurations can be stable minima of our variational 
principle [U Chapter 5.5]. The ad- hoc procedure of working with distributions is justified 
in the paper [3], which also gives concrete hints on how the regularized fermionic projector 
should look like on the Planck scale. For a more detailed stability analysis in the continuum 
see [6]. 
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Despite these results, many questions on the relation between discrete space-time and 
the continuum theory remain open. In particular, it seems an important task to comple- 
ment the picture coming from the discrete side; that is, one should analyze large discrete 
systems and compare the results with the continuum analysis. Since minimizing the action 
for a discrete system can be regarded as a problem of non-linear optimization, numerical 
analysis seems a promising method. Numerical investigations have been carried out suc- 
cessfully for small systems involving few particles and space-time points [7J. However, for 
large systems, the increasing numerical complexity would make it necessary to use more 
sophisticated numerical methods or to work with more powerful computers. Therefore, it 
seems a good idea to begin with simplified systems, which capture essential properties of 
the original system but are easier to handle numerically. In this paper, we shall introduce 
such a simplified system. The method is to employ a spherically symmetric and static 
ansatz for the fermionic projector. This reduces the number of degrees of freedom so 
much that it becomes within reach to simulate systems which are so large that they can 
be compared in a reasonable way to the continuum. 

The paper is organized as follows. In Section [2] we review the mathematical framework 
of the fermionic projector in discrete space-time and introduce our variational principle. 
In Section [3] we take a spherically symmetric and static ansatz in Minkowski space and 
discretize in the time and the radial variable to obtain a two-dimensional lattice. In 
Section [J] our variational principle is adapted to this two-dimensional setting. In Section [5] 
we give a precise definition of our model and discuss its basic properties; for clarity this 
section is self-contained and independent of the rest of the paper. In Section[6]the existence 
of minimizers is proved. In Section[7Jwe present first numerical results and discuss an effect 
of a spontaneous symmetry breaking. We point out that the purpose of this paper is to 
define the model and to discuss some basic properties. Numerical simulations of larger 
systems will be presented in a forthcoming publication. 

2 A Variational Principle in Discrete Space-Time 

We briefly recall the mathematical setting of discrete space-time and the definition of 
our variational principle in the particular case of relevance here (for a more general in- 
troduction see [2]). Let H be a finite-dimensional complex vector space endowed with a 
non-degenerate symmetric sesquilinear form <.|.>. We call (H, <.|.>) an indefinite inner 
product space. The adjoint A* of a linear operator A on H can be defined as in Hilbert 
spaces by the relation <A^ | <!>> = <^ \ A*&>. A selfadjoint and idempotent operator 
is called a projector. To every element x of a finite set M = {1, . . . , m} we associate a 
projector E x . We assume that these projectors are orthogonal and complete, 



Furthermore, we assume that the images E X {H) C H of these projectors are all four- 
dimensional and non-degenerate of signature (2,2). The points x 6 M are called discrete 
space-time points, and the corresponding projectors E x are the space-time projectors. The 
structure (H, <.|.>, (E x ) xG m) is called discrete space-time. Furthermore, we introduce the 
fermionic projector P as a projector on a subspace of H which is negative definite and of di- 
mension /. The vectors in the image of P have the interpretation as the occupied quantum 
states of the system, and / is the number of particles. We refer to (H, <.|.>, (E x ) x£ m, P) 
as a fermion system in discrete space-time. 
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When forming composite expressions in the projectors P and {E x ) x ^Mi it is convenient 
to use the short notations 

= E X V and P(x,y) = E x P E y . (2) 

Using ([1]), we obtain for any <J> G H the formula 

| $> = <*( x ) | $(x)> Ex{H) , (3) 

and thus the vector G E X {H) C P can be thought of as the "localization" of the 

vector $ at the space-time point x. Furthermore, the operator P(x,y) maps E y (H) C H 
to E X (H), and it is often useful to regard it as a mapping only between these subspaces, 

P(x,y) : E y (H) P^P) . 

Again using ([1]), we can write the vector P^ as follows, 

(P*)(x) = P.P* = ^£,P£,* = ^(^P^ftf), 

and thus 

(P*)(x) = ^P(x,y)*(y). (4) 

yeM 

This relation resembles the representation of an operator with an integral kernel. There- 
fore, we call P(x, y) the discrete kernel of the fermionic projector. 

To introduce our variational principle, we define the closed chain A xy by 

A xy = P(x,y)P(y,x) : E X (H) - E X (H) . (5) 

Let Ai, . . . , A4 be the zeros of the characteristic polynomial of A xy , counted with multi- 
plicities. We define the Lagrangian by 

£[4r»] = \ jz^\-\H) 2 (6) 

and introduce the action by summing over the space-time points, 

5[P] = J] . (7) 

x,y€M 

Our variational principle is to minimize this action under variations of the fermionic projec- 
tor. We remark that this is the so-called critical case of the auxiliary variational principle 
as introduced in [H [2] . 

3 The Spherically Symmetric Discretization 

Recall that in discrete space-time, the subspace E X {H) C H associated to a space-time 
point x G M has signature (2,2). In the continuum, this vector space is to be identified 
with an inner product space of the same signature: the space of Dirac spinors at a space- 
time point iEl 4 with the inner product \P<I>, where \P = denotes the adjoint spinor. 
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For any 4x 4-matrix B acting on the spinors, the adjoint with respect to this inner product 
is denoted by B* = 7°S^7°. Furthermore, the indefinite inner product space (H, <-|->) 
in the continuum should correspond to the space of Dirac wave functions in space-time 
with the inner product 

<*|$> = J ^{x)<$>{x)d A x . (8) 

This resembles Q, only the sum has become a space-time integral integral. Likewise, 
in the sum should be replaced by an integral, 



(P*)(x) = J P(x,y)V(y)d 4 y, 



where now P(x, y) is the integral kernel of the fermionic projector of the continuum P. 
Since we assume that our system is isotropic, it follows that it is homogeneous in space. 
Furthermore, we assume that our system is static, and thus the integral kernel depends 
only on the difference y — x, 

P(x, y) = P(f ) for all x, y € M 4 and f := y - x. (9) 

We take the Fourier transform in £, 

/d 4 k ~ 

(2^F P{p) ^ ' (10) 

where p£ denotes the Minkowsi inner product of signature (H ). Let us collect some 

properties of P{p). First, the operator P should be symmetric (= formally self-adjoint) 
with respect to the inner product (JSj) . This means for its integral kernel that 

p(0* = p(-0 , (ii) 

and likewise for its Fourier transform that 

Pip)* = Pip) ■ 

Assuming as in [H §4.1] that the fermionic projector has a vector-scalar structure, P can 
be written as 

P{p) = vjip)^ + 4>(p)l (12) 

with a real vector field v and a real scalar field <j). Moreover, the assumption of spherical 
symmetry implies that the above functions depend only on uj := p° and on k := \p\, and 
that the vector component can be written as 

» i * (1 * k -.i k PI 

VjY = vq 7 + Vk 7 with 7 : = — 

\p\ 

and real- valued functions vq and v^. Next we can exploit that the image of P should 
be negative definite. Moreover, since P should be a projector, it should have positive 
spectrum. Since in Fourier space, P is simply a multiplication operator, we can consider 
the operator Pip) for any fixed p. This gives rise to the conditions that the vector field v 
must have the same Lorentz length as <fr and must be past-directed, 

£'o < and Vq — vj, = (p , 
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and furthermore that <j> must be non-negative. Combining the above conditions, we con- 
clude that P can be written in the form 

P(p) = 4>(w,k) (l - 7° coshr(cj,/c) + 7 fc sinhr(w,A;)) (13) 



with a non- negative function eft and a real function r. Note that we have not yet used 
that P should be idempotent, nor that the rank of P should be equal to the number of 
particles /. Indeed, implementing these conditions requires a more detailed discussion, 
which we postpone until the end of this section. 

We next compute the Fourier transform of (113|> . very similar as in j4j Lemma 5.1]. 
Introducing in position space the polar coordinates £ = (t, r, $, (p) and assuming that r ^ 0, 
the scalar component becomes 



2tt 



<j,(t,r) = ——j I did I k 2 dk dcostf d<p0(u,k)e iujt - ikrcos,& 
( 27r ) J-oo Jo J -i Jo 

oo 

did e iujt I k dk sin(Ax) <f>(id, k) . 
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Vo(t,r) = — —, I duj e lujt I k dk sin(fcr) (f>(id, k) coshr(w,/c) 



-oo 

The zero component of the vector component is computed similarly, 

47T 3 

For the calculation of the radial component, we first need to pull the Dirac matrices out 
of the integrals, 

' **7* = tI^ I A v k (u, k) I 



(2vr) 4 (2vr) 4 7 -i"«v~>"v fe 



1 



oo 



i 7 V ( - ^ J dwe lujt J dk sin(Ax) <f>(id,k) sinh r(u;,/c) 
r jut f^n^ sin(fer) 



47T 3 r ./_oo 'n V kr 



dide luJ / kdk I cos(kr) 4>(id,k) sinhr(a;, k) , 



where we set 7 r = (£t*)/|£|- Combining the above terms, we obtain 

-i roo poo r 

P{g) = — - j due** / kdk<f>(id,k) lsin(fcr) 
47T r J-oo J Q L 



7° coshr(u;,A;) sin(A;r) + i 7 r sinhr(w,/c) (cos(fer) — - — ) . (14) 

V kr J - 



sin(Ar) 
A' ? 

Note that this formula has a well defined limit as r \ 0, and thus we set 

1 poo f'OO 

P(t, r = 0) = j-g / ckde iujt k 2 dk <j>(tjj, k) [l - 7 cosh r(u;, k)] . (15) 

1/ - — 00 J 

In (|13|) and (|14|) . the factors 7 fc and 7 r involve an angular dependence. But all the 
other functions depend only on the position variables (t, r) and the corresponding mo- 
menta (u>, k). We now discretize these variables. In view of (j lip it suffices to consider the 
case t > 0. The position variables should be on a finite lattice £, 

(t,r) G £ := {o,A t ,..., (iV t - l)A f } x {o, A r , . . . , (iV r - 1) A r } , 
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Figure 1: Example of a three-particle system on a 3 x 3- lattice 



where Nt and N r denote the number of lattice points in time and radial directions, 
and At,A r > are the respective lattice spacings. The momentum variables should 
be on the corresponding dual lattice £, 



(u,k) G £ := { - (N t - l)A w ,...,-A w ,o} 



A k ,...,N r A 



(16) 



where we set 



2vr 
A t Nt 



A, 



2vr 



We point out that the parameter uj in (|16|) is non-positive; this is merely a convention 
because we are always free to add to w a multiple of iVjA w . Furthermore, note that 
the points with k = have been excluded in £. This is because the integrands in (I14p 
and (|15p vanish as k \ 0, and thus it seems unnecessary to consider the points with k = 0. 
However, since -P(£) has a non-trivial value at r = (see (I15p ). it seems preferable to take 
into account the points with r = in the lattice £. Replacing the Fourier integrals by a 
discrete Fourier sum, (fT4l) and (fl~5l) become 



p(0 



(1 — 7° coshr) sin( 



(u,fc)e£ 



+ i7 r sinhr (cos(A;r) 



sin(/cr) 
kr 



P(t,r = 0) 



47T 3 



e iaJ 'fc 2 ^ (1-7° coshr), 



ifr/0 (17) 
(18) 



(u,fc)e£ 



with functions and r defined on £. 

The points of the dual lattice £ have the interpretation as the quantum states of 
the system, which may or may not be occupied by fermionic particles. More precisely, 
if 4>(uj,p) 0, a whole "shell" of fermions of energy u and of momenta k with \k\ = p is 
occupied. For most purposes it is convenient and appropriate to count the whole shell of 
fermions as one particle of our lattice model. Thus if (p(uj,p) ^ 0, we say that the lat- 
tice point (u),p) is occupied by a particle; otherwise the lattice point is not occupied. 
A system where n lattice points are occupied is referred to as an n-particle system. 
Each particle is characterized by the values of <j> and r, or, equivalently, by the vector 
(—2k(f> coshr, 2k(f> sinhr). It is convenient to describe the fermion system by drawing these 
vectors at all occupied lattice points, as shown in Figure [T] for a three-particle system. 
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We conclude this section by a discussion of what the parameter / and the idempotence 
condition P 2 = P of discrete space-time mean in the setting of our lattice model. In 
discrete space-time, the number of particles / equals the trace of P. Computing the trace 
of P naively for our lattice model, our homogeneous ansatz ([9]) yields 

/ = TrP = / trP(x, x) d A x = trP(£ = 0) • oo , (19) 

where "tr" denotes the trace of a 4 x 4-matrix. According to (|18p . 

trP(£ = 0) = £ k 2 ^,k), (20) 

(w,fc)e£ 

showing that (|19p is equal to +oo unless P vanishes identically. Here we used essentially 
that, although £ = y — x was discretized on a finite lattice, the space-time variable x 
itself is still an arbitrary vector in Minkowski space. In other words, our lattice system is 
a homogeneous system in infinite volume, and in such a system the number of particles 
is necessarily infinite. The simplest way to bypass this problem is to note that for a 
homogeneous system in discrete space-time [2j Def. 2.4], 

f = ^xP) = ™Tr(PxP), 

and so the number of particles grows linearly with the number of space-time points. Due 
to this simple connection, we can disregard / and consider instead the local trace. This 
has the advantage that the local trace can be identified with the expression (|20|) of our 
lattice system. For the variational principle in discrete space-time dH[7|), it is important 
that variations of P keep the number of particles / fixed. This condition can be carried 
over to our lattice system, giving rise to the so-called trace condition (TC): 

(TC) When varying the fermionic projector of the lattice system (|17[ I18p . the local trace 
as defined by 

/loc := W o k Y] k 2 (f)(uj,k) 

should be kept fixed. 

We conclude that, although / is infinite for our lattice system, the local trace f\ oc is well 
defined and finite. This all we need, because with (TC) we have implemented the condition 
corresponding to the condition in discrete space-time that / should be kept fixed under 
variations of P. We point out that neither / nor /i oc coincides with the number of particles 
as obtained by counting the occupied states. 

The idempotence condition P 2 = P is satisfied if and only if the fermionic wave func- 
tions are properly normalized. As explained above, our lattice model is defined in infinite 
space-time volume, and thus a-priori the normalization integrals diverge. As shown in [TJ 
§2.6], a possible method for removing this divergence is to consider the system in finite 
3-volume and to smear out the mass parameter. However, there are other normalization 
methods, and it is not clear whether they all give rise to the same normalization condi- 
tion for our lattice model. The basic difficulty is related to the fact that each occupied 
lattice point (uj,p) £ £ corresponds to a whole shell of fermions (see above). Thus the 
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corresponding summand in (1171 fl8|) involves an "effective wave function" describing an 
ensemble of fermions. But it is not clear of how many fermions the ensemble consists and 
thus, even if we knew how to normalize each individual fermion, the normalization of the 
effective wave function would still be undetermined. This problem becomes clear if one 
tries to model the same physical system by two lattice models with two different lattice 
spacings. Then in general one must combine several occupied lattice points of the finer 
lattice to one "effective" occupied lattice point of the coarser lattice. As a consequence, the 
normalization of the coarser lattice must be different from that on the finer lattice. This 
explains why there is no simple canonical way to normalize the effective wave functions. 

Our method for avoiding this normalization problem is to choose the normalization in 
such a way that the fermionic projector of the continuum can be carried over easily to 
the lattice system: In Minkowski space, a Dirac sea in the vacuum is described by the 
distribution (see [D §2.2]) 

P(p) = (J + m) 5(p 2 - m 2 ) 0(-p°) . (21) 

Taking the Fourier transform and carrying out the angular integrals, we obtain again the 
expressions (|14[ I15p . but now with <p{oj,k) = 5(lo 2 — k 2 — m 2 ). This allows us to carry out 
the fe-integral, 



oo roc 



doj e iujt I k dk 5(co 2 - k 2 - m 2 ) ■ ■ ■ = \ du) e iujt | 



oo JO 



k=yu 



The easiest method to discretize the obtained expression is to replace the w-integral by a 
sum, and to choose for every u < —m a lattice point (lj, k) G £ such that 



< k - \]u 2 - m 2 < A k . (22) 

An example for the resulting discretized Dirac sea is shown in Figure El 

Note that for this configuration, <j){u),k) = l/(2k) at all occupied lattice points. Next 
we allow to modify this configuration, as long as the normalization integrals remain un- 
changed: First, changing r corresponds to a unitary transformation of the corresponding 
state, without influence on the normalization. Second, hopping from a lattice point to 
another unoccupied lattice point with the same value of k changes the state only by the 
phase factor exp(— i{uj — uj')t), again without influence on the normalization. This leaves 
us with the so-called normalization condition (NC): 
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(NC) The function <\> in (fTTl [TH]) should only take the two values 

4>(w,k) = or <t>{u>,k) = — . 

We again point out that this normalization condition is not canonical. It could be modified 
or even be left out completely. It seems an interesting question to analyze how the behavior 
of the lattice model depends on the choice of the normalization condition. 



4 The Variational Principle on the Lattice 

The Lagrangian ([6]) is also well defined for our lattice model. Let us compute it in more 
detail. We decompose the fermionic projector (|144 [15]) into its scalar and vector compo- 
nents, 



i>(0 = <t>{t,r) l + v (t,r)~/° +v r (t,r)Y = 01 + 



VjJ 



j 



Furthermore, using that the functions <f> and r are real, we find that 

P(-0 = P(0* =^l + vj7 j - 
Thus, omitting the argument £, the closed chain © becomes 

A = y + </>)$ + ■$). 

For the computation of the spectrum, it is useful to decompose A in the form 

A = A 2 + A 2 + n 

with 



A, = i 

2 



A2 = 4> i + i> 4> , n = Vjvi + 4>4> . 

A short calculation shows that the matrices A\ and A2 anti-commute, and thus 

(A-fi) 2 = A\ + A 2 2 = D[A] 1 , (23) 

where we set 

D[A] = - tr(A 2 ) - — (trA) 2 = {v j v3) 2 -\v j v j \ 2 + {v j + (pv]){v j ^ + (f> : rf). (24) 



The identity (|23|) shows that the characteristic polynomial of the matrix A has the two 
zeros 

A ± = Vj vi+(f>0±VD. (25) 

If these two zeros are distinct, they both have multiplicity two. If the two zeros coincide, 
there is only one zero of multiplicity four. Hence the Lagrangian ([6]) simplifies to 

C[A] = (|A + |-|A_|) 2 (26) 

In order to further simplify the Lagrangian, we introduce a discrete causal structure, in 
agreement with [5]. 
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Definition 4.1 A lattice point (t,r) € £ is called 



timelike if D[A(t,r)] > 
spacelike if D[A(t, r)] < 0. 

If (t, r) is spacelike, the X± form a complex conjugate pair, and the Lagrangian 
vanishes. If conversely the discriminant is non-negative, the \± are both real. In this 
case, the calculation 

A + A_ = (vv + 44f - [{vvf -v 2 v 2 + {v~4> + 4>vf] 
= 2{vv)\4>\ 2 + \4>\ i + v 2 v 2 - (v$ + <f>v) 2 

= \^ + v ^y2_ v 2^_ (j)2 _ 2 = {v 2_^2 ){ -2_^ ) > q 

(where we omitted the tensor indices in an obvious way) shows that A_|_ and A_ have the 
same sign, and so we can leave out the absolute values in (I26p . We conclude that 

, , ( AD[A(t,r)] if (t, r) is timelike 
1 otherwise , 

where D is given by (|24|) , Hence our Lagrangian is compatible with the discrete causal 
structure in the sense that it vanishes if (t, r) is spacelike. 

Before we can set up the variational principle, we need to think about what the sum 
over the space-time points in ([7]) should correspond to in our lattice system. Since we are 
considering a homogeneous system, one of the sums simply gives a factor m, and we can 
leave out this sum. The other sum in the continuum should correspond to a space-time 
integral (see for example ©)• In our lattice system, the point (t,r) can be thought of as 
the 2-dimensional sphere |£| = r at time t. Therefore, we replace the spatial integral by 
a sum over the discretized radii, but with a weight factor which takes into account that 
the surfaces of the spheres grow quadratically in r. More precisely, we identify (t, r) with 
a shell of radius between r — A r /2 and r + A r /2. This leads us to the replacement rule 

Nr-l 



f d? ... A 3 r V Pr (nA 



n=Q 

with the weight function p r given by 



4vr J 1/8 ifn = 

" ~3~" I (n + 1/2) 3 -(n-l/2f if n > . 



When discretizing the time integral, we need to take into account that on the lattice £, 
the time parameter t is always non-negative. Since the Lagrangian is symmetric, C[A xy ] = 
£[A yx ] (see [H §3.5]), this can be done simply by counting the lattice points with t > 
twice. Thus we discretize the time integral by 



/oo 
-oo 



JVr-1 



dt ... — ► A t ^ Pt(nA r ) 

J — oc 

with 



n=0 
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Then the action becomes 



S[P] = A t A T 3 Pt(t)Pr(r)C(t,r). 

(t,r)e£ 

Our variational principle is to minimize this action by varying the functions <fi and r in (|17l 
TSj) under the constraints (TC) and (NC). 

With the constructions of Sections [3] and H] we successively derived our two-dimensional 
lattice model. Clearly, not all the arguments leading to the model were rigorous, and also 
we put in strong assumptions on the physical situation which we have in mind. More pre- 
cisely, the main assumption was the spherically symmetric and static ansatz with a vector- 
scalar structure (jlOl [T2]) ; this ansatz was merely a matter of convenience and simplicity. 
Moreover, the choice of the weight function p involved some arbitrariness. However, we 
do not consider this to be critical because choosing the weight factors in (|27p differently 
should not change the qualitative behavior of the model (except that for the existence 
of minimizers it is important that p r (0) ^ 0; see Section [6]). Finally, the normalization 
condition (NC) could be modified, as discussed in detail at the end of Section El 

The main point of interest of our lattice model is that it allows to describe discretiza- 
tions of Dirac seas f)22f) but also completely different configurations of the fermions. Thus 
within the lattice model it should be possible to analyze in detail whether and how Dirac 
sea configurations form as minimizers of our variational principle. Moreover, in our lat- 
tice model one can implement all the spherically symmetric regularization effects as found 
in Hence our lattice model should make it possible to verify effects from [3] coming 
from the discrete side and to analyze these effects in greater detail. 

In the next section we shall define our lattice model once again more systematically, 
making the following simplifications: 

• By scaling we can always arrange that A w and A^ have an arbitrary value. It is 
most convenient to choose 



A* = 1, A fc = 1. 

Then 

A t = — , A. r = — 

N t N r 



The formulas for P, (\17\ I18 p . only involve the two Dirac matrices 7 and 7 r , which 
satisfy the anti-commutation rules 

(7°) 2 = 1, iff = -1, {7°,7 r } = 0. 

Since these anti-commutation rules can be realized already by 2 x 2-matrices, we 
may simplify the matrix structure by the replacements 

7° — a 3 , 7 r — 

where a 1 are the usual Pauli matrices. Modifying the definition of the discrimi- 
nant ([MD to 

D[A] = Itr^ 2 )-^) 2 

(where now "tr" clearly denotes the trace of a 2 x 2-matrix) , the Lagrangian remains 
unchanged. 
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• In order to simplify the normalization condition (NC), it is convenient to introduce 
the function 

$(k,u) = 2k(f)(k,Lo) , 
where for notational simplicity we also omitted the tilde. 

• In order to simplify the prefactors, we multiply P by 4tt, divide the Lagrangian by 
four, and divide the action by (27r) 4 7r/6. Furthermore, we multiply f\ oc by a factor 

Of 7T 3 . 



5 Definition of the Model and Basic Properties 

For given integer parameters N t , N r and p we introduce the lattice £ and its dual lattice £, 

. . r 2vr N t - 1 1 r 2vr N r - 1 1 

ftr) 6 fi = {0, 2,—} x{0,-....,2»_} 

(u,k) € £ = {-(JV,-1), -l,o} x {l, JV r } . 

On £ we choose a non- negative function <E> and a real function r, which vanish except at p 
lattice points. We set 



P(f) = - 5^ e icj< $ [(1-a 3 coshr) sin(£;r) 



r 

1 ( s sin(A;r)\i 

+ 0- 1 sinhr cos(Ax r — " , if r 7^ (29 

V kr / J 

P(i, r = 0) = ^ e^* k $ (1 - a 3 coshr) , (30) 

(u;,fc)e& 

where a 1 and cr 3 are two Pauli matrices. For any (t, r) 6 £ we introduce the closed 
chain x4(t, r) by 

A(t,r) = P(t,r)P(t,r)* , 
where the adjoint with respect to the spin scalar product is given by 

P(t,r)* = a 3 P(t, r y a 3 , 

and the dagger denotes transposition and complex conjugation. We define the discrimi- 
nant D[A] and the Lagrangian C[A] by 

D[A] = \ tr(x4 2 ) - \ {t?Af (31) 
C[A] = D[A] @(D[A]) , (32) 

where is the Heaviside function. The action is 

5 = Njp E Pt(t) Pr (r)C[A(t,r)], 

1 r (t,r)6fl 
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where pt and p r are the weight functions 

n 

n_\ i 1 ifn = 

27r W r ) ~ \ (2n + l) 3 -(2n-l) 3 if n > . 

Our variational principle is to minimize the action, varying the functions <1> and r under 
the following constraints: 

(TC) The local trace 

/ loc := £ k$(u,k) (33) 

should be kept fixed. 

(NC) The function <1> should only take the two values &(u>, k) = or 3>(cj, A;) = 1. 

The last condition (NC) could be weakened or left out (see the discussion at the end of 
Section Ej). 

According to Definition 14.11 the functions <J> and r induce on £ a discrete causal 
structure. The Lagrangian is compatible with the discrete causal structure in the sense 
that it vanishes if (t, r) is spacelike. Furthermore, our lattice system has the following 
symmetries: 

symmetry under parity transformations: The traces in f)31 1) vanish unless an even 
number of matrices a 1 appears. Therefore, the Lagrangian remains unchanged if the 
factor sinh r in (|29p flips sign. Hence the action is symmetric under the transforma- 
tion 

r(a/,fc) — ► -r(u,k) for all (u, k) G £. (34) 

This transformation changes the sign of the spatial component of P. The name 
"parity transformation" comes from the analogy to the usual parity transforma- 
tion x — ► — x. 

gauge symmetry: We introduce on the dual lattice £ for any !! £ Z the translation 
respecting the periodic boundary conditions 

u — ► Co = (oj + fi) mod JVt (35) 

and also translate the functions r and $ by setting 

f(£b,k) = t(u, k) , $(Cj,k) = $(u),k) . 

This translation in momentum space corresponds to a multiplication by a phase 
factor in position space, 

P(£) = e J ™ P(£) . 

This phase factor drops out when forming the closed chain, and thus the Lagrangian 
remains unchanged. The transformation (|35p are precisely those local gauge trans- 
formations which are compatible with our spherically symmetric and static ansatz. 
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6 Existence of Minimizers 



In this section we prove an existence result, which is so general that it applies also in the 
case when the normalization condition (NC) is weakened. 

Proposition 6.1 Consider the variational principle of Section [5| with the trace condi- 
tion (TC) and, instead of (NC), the weaker condition that that there is a parameter e > 
such that 

$(u>,k) = or $(u,k) > e for all (u , k) £ £ . 
Then the minimum of the action is attained. 

Proof. Since the Lagrangian is non-negative, we can estimate the action from above by 
the Lagrangian at the origin t = r = 0, 

S > £[A(0,0)] . (36) 

At the origin, the fermionic projector takes the form (see (|30p ) 

P(0) = £;$(l-o- 3 coshr). 

(w,fc)e£ 

This matrix can be diagonalized and has the two eigenvalues 

(i± = ^2 k $ (1 ± coshr) . 

Thus the closed chain A(0,0) has the two eigenvalues A± = As a consequence, 

using d3U) and US}, 

= \ + (i-) 2 (M+-M- ) 2 
\ 2 

(37) 

Consider a minimal sequence. Then, according to (|36p . the expression (|37p is uniformly 
bounded. If /i oc = 0, our system is trivial, and thus we may assume that /i oc is a positive 
constant. Using (|33p and the fact that k > 1, we conclude that the functions $ are 
uniformly bounded. The boundedness of (|37p implies that there is a constant C > such 
that 

A; $ cosh t < C 

(w,fe)6£ 

for all elements of the minimal sequence. Whenever $ vanishes, we can also set r equal 
to zero. If <3? is non-zero, the inequality $ > e gives a uniform upper bound for coshr, 

coshr < — . 

e 

We conclude that the functions $ and r are uniformly bounded. Hence a compactness 
argument allows us to choose a convergent subsequence. Since our action is obviously 
continuous, the limit is the desired minimizer. ■ 

We point out that this proposition makes no statement on uniqueness. There seems no 
reason why the minimizers should be unique. In Section [7] we shall see examples with 
several minimizers. 
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Figure 3: Action for a two-particle system on an 8 x 6-lattice 

7 First Numerical Results, A Mechanism of Spontaneous 
Symmetry Breaking 

We now discuss numerical results for a simple lattice system, with the intention of ex- 
emplifying a few general properties of our model. More precisely, we choose Nt = 8, 
N r = 6 and occupy two lattice points, one with momentum k = 1 and one with k = 2. 
Satisfying the condition (TC) and (NC), the system is characterized by the two discrete 
parameters Ui,U2 £ {— 7, ...,0} and the two real parameters n and T2 at the occu- 
pied lattice points. We find that the absolute minimum of the action is attained if we 
choose ui = — 1 and u>2 = —2. For this choice of the occupied lattice points, Figure [3] 
shows a contour plot of the action in the (t±, T2)-plane. The plot is symmetric under 
reflections at the origin; this is the symmetry under parity transformations (|34H . The 
minimum at the origin corresponds to the trivial configuration n = T2 = 0, where the 
two vectors v- t := (— $j coshrj, $j sinhrj) are both parallel to the w-axis. However, this 
is only a local minimum, whereas the absolute minimum of the action is attained at the 
two points (t\ ~ 1.5, r% ~ 1) and (n ~ — 1.5, ~ — 1). The interesting point is that 
the absolute minima are non-trivial, meaning that the occupied points distinguish specific 
lattice points and that the corresponding vectors Vi are not parallel. The minimizers are 
not symmetric under (|34p . Thus by choosing one of the minimizers, the symmetry under 
parity transformations is spontaneously broken. 

This effect of spontaneous symmetry breaking can be understood in analogy to the 
Higgs mechanism in the standard model, where the double-well potential has non-trivial 
minima, and by choosing such a minimizer the corresponding vacuum breaks the original 
SU (2)-symmetry. For our lattice system, the role of the Higgs potential is played by the 
particular form of our action in combination with the constraints as given by the trace 
condition and the normalization condition. 

We remark for clarity that the effect of spontaneous symmetry breaking observed here is 
much different from the effect described in [3] for general fermion systems in discrete space- 
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time. Apart from the fact that in [3] we consider instead of parity the outer symmetry 
group, the main difference is that in [3] we do not specify the action, but the spontaneous 
symmetry breaking arises merely by constructing fermionic projectors for a given outer 
symmetry group. Here, on the contrary, the specific form of the action is crucial. This 
is made clear by the alternative action S = k)e£ \ T ( UJ i ^)l> wmcn only has trivial 
minimizers with r = which are all symmetric under the parity transformation. 

We conjecture that for a large lattice and many particles there are minimizers which 
look similar to the discretized Dirac sea structure of Figure [2j Since the presentation of 
larger simulations is more involved and also requires a detailed description of the numerical 
methods, we postpone this analysis to a forthcoming publication. 
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